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ABSTRACT 

Reconnection of alternating magnetic fields is an important energy transfor- 
mation mechanism in Poynting dominated outflows. We show that the recon- 
nection is facilitated by the Kruskal-Schwarzschild instability of current sheets 
separating the oppositely directed fields. This instability, which is a magnetic 
counterpart of the Ray leigh- Taylor instability, develops if the flow is accelerated. 
Then the plasma drips out of the current sheet providing conditions for rapid 
reconnection. Since the magnetic dissipation leads to the flow acceleration, the 
process is self-sustaining. In pulsar winds, this process could barely compete 
with the earlier proposed dissipation mechanisms. However, the novel mecha- 
nism turns out to be very efficient at AGN and GRB conditions. 

Subject headings: galaxies:jets - gamma-rays:bursts - MHD - relativity 

1. Introduction 

It is now widely agreed that relativistic outfiows are launched hydromagnetically. The 
commonly invoked scenario is that a strong magnetic field in a rapidly rotating compact 
object serves to convert the rotational energy into the energy of the outfiow (see, e.g., the 
review by Spruit (2010)). The main advantage of the magnetic launch mechanism is that the 
magnetic field lines, like driving belts, could transfer the rotational energy to a low density 
periphery of the central engine thus forming a baryon pure outfiow with the energy per 
particle significantly exceeding the rest mass energy. Such outfiows potentially capable of 
reaching relativistic velocities. Since the energy is transported, at least initially, in the form 
of the Poynting fiux, the question arises of how and where the electromagnetic energy is 
eventually transformed to the plasma energy. An important point is that the fiow could be 
considered as true matter dominated only when the electro-magnetic energy falls well below 
the equipartition level because only then the efficient shock dissipation becomes possible 
(Kennel & Coroniti 1984). 
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In the scope of ideal magnetohydrodynamics (MHD), the energy could in principle be 
transferred to the plasma via gradual acceleration by electromagnetic stresses. The problem 
is that in rclativistic flows, the magnetic and electric forces nearly compensate each other 
so that the flow is nearly ballistic. It turns out that unconfincd flows stops accelerating 
when still being highly Poynting dominated. A more efficient acceleration regime is achieved 
if the flow is coUimated by an external medium. Such a configuration arises naturally in 
gamma-ray bursts (GRBs), where the relativistic jet from the collapsing stellar core pushes 
its way through the stellar envelope. In accreting systems, the outflow from the rotating 
black hole could be coUimated by the pressure of a slow (and generally magnetized) wind 
from the outer parts of the accretion disc. However, a systematic study shows (Komissarov 
et al. 2007, 2009; Tchekhovskoy et al. 2008, 2010; Lyubarsky 2009) that even though 
externally conflned flows could in principle be accelerated till the equipartition between the 
kinetic and magnetic energy, the conditions of that are rather restrictive. Moreover, the true 
matter dominated stage could in any case be achieved only at exponentially large distances 
(Lyubarsky 2010). 

An efficient conversion of the electro-magnetic into the plasma energy was claimed for 
impulsive flows (Granot et al. 2010,; Lyutikov 2010a, b). In a freely expanding shell, the 
magnetic energy is indeed eventually transformed into the kinetic energy. However, Levinson 
(2010) demonstrated that in the real world, it is difficult to provide conditions for the free 
expansion because even a tenuous ambient medium hampers the expansion. Intermittent 
ejection of small sub-shells could hardly help because the shells merge while still highly 
magnetized unless the distance between them is much larger than their width. 

The problem of the energy conversion could be resolved if one abandons the assump- 
tion of no dissipation. However, one has to find a reliable dissipation mechanism in highly 
conductive space plasma. One of the ideas is that some sort of MHD instabihty destroys 
the regular structure of the magnetic field thus triggering the anomalous dissipation. The 
kink instability is the best candidate (Begelman 1998; Giannios & Spruit 2006). However, 
any global MHD instability could develop only if the proper Alfven crossing time is less than 
the propagation time, which implies 9^ < 1, where 7 is the jet Lorentz factor, 6 the opening 
angle. This condition is rather restrictive; for example it could hardly be fulfilled in GRBs 
(e.g. Tchekhovskoy et al. 2010). Even if the kink instability does develop, it is still not clear 
whether the flow is disrupted or just helically distorted. The impact of the instability on 
the jet structure could be studied only with three-dimensional numerical simulations. Till 
now only a few attempts was made, the results still being controversial (Moll & Spruit 2008; 
(McKinney & Blandford 2008). 

Magnetic dissipation could occur if a small-scale structure with oppositely directed 
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fields preexisted in the flow. In pulsar winds, such a structure arises naturally because the 
magnetic field in the equatorial belt of the flow changes sign every half of period. It is 
the annihilation of the oppositely directed fields that provides the main energy conversion 
mechanisms in pulsar winds (Coroniti 1990; Lyubarsky & Kirk 2001; Kirk & Skj^raasen 
2003; Petri & Lyubarsky 2007; Zenitani & Hoshino 2007; see also the review by Kirk et al. 
(2009)). The outflow with alternating fields could in principle arise also in accreting systems 
if the magnetic field in the central engine changes sign. Then independently of the exact 
field configuration at the launch site, the flow expansion ensures that in the far zone, the 
overall magnetic structure is that of the "striped wind" with stripes of oppositely directed 
azimuthal fields separated by current sheets. The width of the stripes is small as compared 
with the scale of the flow therefore locally the structure of the flow is very simple: plane 
current sheets separate domains with oppositely directed magnetic fields. Within the current 
sheets, the magnetic field changes sign and the pressure balance is maintained by the thermal 
pressure of the plasma. 

Efficient annihilation of oppositely directed magnetic fields across the current sheet 
implies some sort of strong anomalous resistivity. The necessary condition is that the drift 
velocity of the current carriers becomes large enough, which is equivalent to the condition 
that the particle Larmor radius is not small as compared with the thickness of the sheet. 
Since the current in the sheets decays as whereas the plasma density drops down as 
r~^, the anomalous dissipation due to the charge starvation should inevitably occur at a 
large enough distance from the source. In pulsar winds, the plasma density is extremely 
low so that the alternating fields do dissipate in the far zone of the wind or at the wind 
termination shock (Kirk & Skjaeraasen 2003; Petri & Lyubarsky 2007; Zenitani & Hoshino 
2007). However in accreting systems, the outfiow is expected to be heavily loaded by the 
plasma from the accretion disk (and in GRBs, also by the plasma from the progenitor star 
(Levinson & Eichlcr 2003)) therefore the charge starvation conditions are achieved only at 
unreasonably large distances. 

Spruit et al. (2001) postulated that in GRB outfiows, the alternating fields annihilate 
with a rate ~ 0.1 of the speed of light. The model based on this assumption is capable of 
explaining many essential features of GRBs (Drenkhahn 2002; Drenkhahn & Spruit 2002; 
Giannios 2006, 2008; Giannios & Spruit 2007). The problem is that the reconnection rate of 
the order of 10% from the Alfven velocity is the maximal possible reconnection rate, which is 
achieved only at special conditions (see review by Yamada et al. 2010). For example, there 
is practically no reconnection across the heliospheric current sheet, which is a prototype of 
large-scale current sheets in any astrophysical outflow. 

In this Letter, we propose a mechanism for fast reconnection of the alternating flelds in 
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Poynting dominated outflows. The basic idea is the following. If the flow is accelerated (or 
decelerated), an effective gravity force appears in the proper plasma frame so that the hot 
plasma within the current sheet is supported against the gravity force by the magnetic pres- 
sure. Such a configuration is known to be unstable with respect to the Kruskal-Schwarzschild 
instability, which is a magnetic counterpart of the Rayleigh- Taylor instability. Under the in- 
fluence of the effective gravity force, the plasma could drip down between the magnetic field 
lines (e.g. Infeld (1989)) so that bridges between trickles thin out until the current sheet tears 
(Fig. 1). Thus the instability facilitates the reconnection. An important point is that the 
magnetic dissipation forces the flow to accelerate (Lyubarsky & Kirk 2001; Drenkhahn 2002; 
Drenkhahn & Spruit 2002; Kirk & Skjseraasen 2003) therefore the process is self-sustaining: 
the acceleration is an aid to the reconnection whereas the reconnection promotes the accel- 
eration. In this Letter, we show that the alternating fields could be completely dissipated 
by this mechanism and find the characteristic dissipation scale. 



2. Kruskal-Schwarzschild instability of the current sheet 

Let us consider a hot, nonmagnetized plasma slab, < z < A, separating domains with 
oppositely directed magnetic fields. Let the field be directed along the y axis. We assume 
that this structure is accelerated in z direction; then an effective downward gravity force 
arises in the proper frame of the flow. In the flow frame, the continuity and Euler equations 
for slow motions are written as 

-- + V-nv = 0; (1) 
h— = -c^Vp + hg; (2) 



where g is an effective gravity. 



h = nine? + — — -p (3) 



the specific enthalpy, F the adiabatic index. 

In the strongly magnetized domain, the flux freezing condition ties the plasma and 
the magnetic field together. Therefore when considering flat motions perpendicular to the 
magnetic field {vy = 0), one can take B/n = const. This implies the effective equation of 
state p = B'^/Stc oc n^; h = B'^/Att = 2p so that one can describe the magnetized domain by 
Eqs. ^ and © with m = and F = 2. 



We consider evolution of small perturbations, h = Hq + hi, p = po + pi, n = + rii, 
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where the unperturbed state satisfies the equihbrium equation 

c — = -gh,. (4) 
We assume that the initial state is nearly homogeneous, which implies 

ghoA < c%. (5) 

Both in the magnetized domain and in the relativistically hot slab, p ^ h therefore this 
condition could be written as 

A « c^g. (6) 

Stability analysis implies harmonic perturbations of the form exp{—iujt + ikx). In this 
case, Eqs. ([1]) and ([2]) are reduced to 

— iurii + no ( ikvx + —7—] = 0; (7) 



dz ^ 

uhov^ = kc^pi, (8) 

-iuhoVz = -c^-r- - gh- (9) 
az 



Eliminating velocities, one gets 



d'^Pi g dhi 2 



dz'^ & dz c^rio 

One can express pi via ni making use of Eq. and the thermodynamic identity d{h/n) = 
dp/n; this yields 

= (r - 1) ( — -mcM ni. (11) 
Now one can write the closed equation for rii. 

Tpo(Pni^ f ^ Po \ drii ^ fuj^hp j^2'^Po\ ^ 

no dz"^ \ r — 1 c'^uq J dz \ no uq J 

Note that we implicitly used the condition (|6]) that the system is nearly homogeneous. 

Looking for the solution in the form rii oc exp kz, one finds 



g f T ^ mnoC^\ ^ g^ / T ^ mnoC^ V ^ hp u"^ ^^^^ 



2c2 Vr - 1 Tpo J V 4c4 Vr - 1 Tpo J Tpo 6 

The condition ([6]) suggests that for kA > 1, one can neglect the terms with g as compared 
with those with k. Qualitative considerations, which will be justified by the exact solution. 
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show that the growth rate is determined by the free-fall time-scale, u ~ \^kg <^ k, therefore 
one can also neglect the term with u. Then Eq. ( fT3l) is reduced simply to n = ±k. 

Taking into account that the perturbation should decay far enough from the slab, one 
can write 

= Cie'''; z<0; (14) 

= C2e^^ + C^e-^'] < ^ < A; (15) 

p^ = de-^^ z>A; (16) 

where the constants C could be found from the condition of continuity at the boundaries of 
the slab. 

Let us consider the boundary z = 0. Defining the plasma shift, ^, as d^/dt = Vz] i.e. 

e = A (17) 
(jj 

one can write the continuity of the flow and of the pressure as 

vM = <m (18) 

Po(0+Pi(o)=Po(0+K(o); (19) 

where the primed quantities are referred to the hot slab and the unprimed ones to the highly 
magnetized medium. For the wavelength kA > 1, Eq. (Q is reduced, with account of the 
condition (jG]), to 

1 dpi 



^ uho dz 

Making use of Eq. ( !T4l) for pi and Eq. (ITSl) for p[, one writes Eq. (fTSjl as 



(20) 



^Ci + ^{C2-C,) = 0. (21) 
ho 

The pressure continuity condition ( IT9|) is reduced to an equation for Ci, C2 and C3 if one 
takes into account that according to Eqs. @, (ITTl) and (l20l) 

Po(0 = MO) + = MO) - -i-,^. (22) 
az c'^oj'^ dz 



Then one gets 



Ci-C2-C^ = 0. (23) 
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Applying the same considerations to the boundary 2; = A, one finds two more equations 

J_ (c^e^^ - Cse-''^) + = 0; (24) 
/^o ho 

C,e'^ + Cae-'^^ - C, (^1 - ^-^Z^^ = 0. (25) 

Now we have four hnear homogeneous equations for the coefficients C. The solvabihty 
condition for this set of equations provides the dispersion equation. After simple algebra, 
one gets 

uj' = ' - ^"P(-2feA) ^26) 

-exp(-2A:A) 

One sees that this equation always possesses an unstable root, Imu < 0. One can check 
that if the gravity force is directed downwards, g > 0, the bottom boundary of the slab is 
unstable whereas g < implies instability of the upper boundary. 

In the strongly magnetized medium, ho = 2po. In the relativistically hot slab, one can 
write /iq = Ap^ = 4po- Then one gets 



ho + hp 
ho - h'r. 



(27) 



In the small wavelength limit, kA ^ 1, the growth rate could be found from equations (1261) 
and (EZD as 

V = -lrnco=^, (28) 

which coincides with the growth rate for the instability of the boundary between the two 
half-spaces, one filled with the plasma and another with the magnetic field. In the opposite 
limit, kA ^ 1, one finds 

r]=(^y^\'/^e/\ (29) 



^2) 

In this regime, the plasma slab is bent as a whole. 



3. Interplay between magnetic dissipation and acceleration 

Now let as consider a highly relativistic outfiow in which the magnetic field forms stripes 
of the opposite magnetic polarity containing cold plasma separated by currents sheets con- 
taining hot plasma. Let us denote the width of the stripe I. In pulsar winds, I is of the 
order of the light cylinder radius; specifically in the equatorial plane, / = vrr^. In accreting 
systems, / could be as small as the gravitational radius, I ^ Vg. 
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When the flow accelerates, the effective gravity in the proper frame is 



2 0^7 



In this section, all quantities measured in the proper plasma frame are marked by primes. 
As it was demonstrated in the previous section, the hot plasma drips out of the sheet in this 
case so that the slab eventually tears at some points (see Fig. 1) facilitating the magnetic 
reconnection. The most dangerous are perturbations with the wavelength comparable with 
the thickness of the slab. Small wavelength perturbations just produce ripples on the slab 
surface whereas long wavelength perturbations result only in bending of the slab. Therefore 
the characteristic reconnection time is determined by the instability growth rate for the 
wavelength of the order of the sheet width. 



r' ~ ^//\'/g'. (31) 

As the field annihilates, the current sheet expands consuming the energy and the plasma 
from the magnetized domain. The rate of expansion could be written as 

dA' A' 



C- = CvW. (32) 

where C is a numerical factor. Three-dimensional numerical simulations of the current sheet 
supported against the gravity force are necessary in order to confirm this scaling and to 
provide an estimate for C^. 

Substituting equation ( l30l) into equation ( l32l) and transforming to the lab frame, r = 
ct'7, A = A'/7, one can write the rate of the sheet expansion as 



^ = CV7A§I. (33) 

dr V dr 

An important point is that dissipation of the Poynting flux forces the flow to accelerate 
(Lyubarsky & Kirk 2001; Drenkhahn 2002; Drenkhahn & Spruit 2002; Kirk & Skj^raasen 
2003) therefore the process of dissipation due to the Kruskal-Schwarzschild instability is 
self-sustaining: the acceleration facilitates the dissipation whereas the dissipation leads to 
the acceleration. At the initial stage, Poynting dominated flows are accelerated even with- 
out dissipation therefore the process could be initiahzed naturally. Below we present the 
self-consistent description of the acceleration/dissipation process making use of the short- 
wavelength approximation, / <^ r, developed by Lyubarsky & Kirk (2001). 

We consider the radial flow with a purely toroidal magnetic field. This approximation 
is valid far enough from the source where one can neglect the poloidal field in the proper 
frame of the flow, 

7 > BjBp = R/tl. (34) 



Here R is the cylindrical radius of the flow. 



The flow is steady on the average therefore the continuity and the energy conservation 
equations could be written as 



r {nv) 



h'j V + - — V 

4:71 



const; 



const: 



(35) 
(36) 



where angle brackets denote averaging over the striped structure. Adopting for simplicity 
the rectangular stripe shape such that within the current sheet, the magnetic fleld is zero 
and the plasma is relativistically hot with the pressure balancing the magnetic pressure of 
the cold, magnetized domain, h' = Ap = B''^/2it, one writes the conservation equations as 



2 2 

f 7 r 



vr'^lucA + {I — A)n/i] = const; 
-A + (/-A) (m<c^ + ^j 



const. 



(37) 
(38) 



Here Uc and are the plasma density in the cold magnetic domain and in the hot current 
sheet, respectively. The constants are determined by the parameters at the inlet of the flow, 
r = ro, which should be chosen far enough from the central source so that the condition ( IMl) 
is already fulflUed. 

In the limit A <^ /, 7 ^ 1, these equations may be combined into the asymptotic 
relation (equation (29) in Lyubarsky & Kirk 2001; note the change in notations: Lyubarsky 
& Kirk used dimensionless A, which is equal to our A/Z) 



Ur 



A 1 
T + 2^ 



1 



7o 



- 1; 



(39) 



where 70 is the Lorentz factor at the inlet of the flow 

fi2 



Anmnc-yc^ 



(40) 



the initial magnetization parameter. This equation should be complemented by the entropy 
equation and by the prescription for the reconnection rate. In our case, the last is given by 
equation ( l33l) . The entropy equation implies Uh = Suc at large distances (Lyubarsky & Kirk 
2001). 



When A/ 1 ^ l/7o) equation (l39l) is reduced to 
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Substituting this relation into equation fl33|) . one gets the equation for 7 



2 (i7 3C Vo7o 



^ (42) 

which is immediately solved to yield 

where 7max = o"o7o is the Lorentz factor achieved if the Poynting flux is completely trans- 
formed into the plasma kinetic energy. Combining this relation with equation (HT]) one sees 
that the Poynting flux is dissipated completely, A ~ /, at the scale 

rdiss = 12(^y/. (44) 

At this scale, 7 approaches 7max- 



4. Discussion 

In this Letter, we proposed a novel mechanism for the energy release in Poynting dom- 
inated outflows. The mechanism operates if a significant fraction of the Poynting flux is 
transferred by alternating fields. Such a structure inevitably arises in the wind from obliquely 
rotating pulsars and magnetars. In accreting sources, one has to assume that the magnetic 
polarity changes sign near the black holes. Earlier only anomalous resistivity mechanisms for 
the field dissipation were considered; they are efficient only in the current starvation regime 
when the thickness of the current sheet is comparable with the Larmor radius of the particles 
(Coroniti 1990; Lyubarsky & Kirk 2001; Kirk & Skjaeraasen 2003; Zenitani & Hoshino 2007). 
We found a universal mechanism for dissipation of alternating fields, which works even if 
the conditions of ideal MHD are fulfilled and anomalous resistivity does not appear. Thus 
we justified the assumption of Spruit et al. (2001) that in Poynting dominated outfiows, the 
annihilation rate of alternating fields comprises a good fraction of the speed of light. 

The mechanism is based on the MHD instability of the current sheet in the accelerating 
fiow. Under the infiuence of the inertia force, the plasma drips out of the current sheet 
allowing the oppositely directed fields come together. The ensuing field annihilation leads to 
the fiow acceleration so that the process is self-sustaining. Note that this instability is easily 
suppressed by magnetic shear therefore it is very important that in the far zone of the fiow, 
the field becomes predominantly toroidal, i.e. practically shearless. 
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We have shown that the complete dissipation occurs at the scale (1441) . In pulsar winds, 
7max ~ 10^ — 10^ and / ~ r^, so that with this mechanism, the alternating fields are dissipated 
at the distance ~ 10^ — lO^^r^, which is comparable or even larger than the dissipation scale 
due to the current starvation mechanisms (Lyubarsky & Kirk 2001; Kirk & Skj^raasen 
2003; Zenitani & Hoshino 2007). However, this mechanism could efficiently work in GRBs 
and AGNs because 7max in these systems is not very large whereas the current starvation 
distance is huge. 

In AGNs, 

7max ~ 10 and I ^ fg therefore this mechanism provides the energy release 
at the scale of ~ lOOOr^. This is compatible with observations because fitting of the blazar 
spectra implies that jets are already matter dominated at ~ lOOOr^ (Ghisellini et al. 2010). 
The proposed mechanism provides also a natural basis for the "jets- in- a- jet" model (Giannios 
et al. 2009, 2010) that accounts for the observed ultra-fast TeV variability of blazars. 

In GRBs, 7max ~ 1000 so that the energy release occurs at the distance of the order of 
lO^'^ cm. The magnetic dissipation at this scale could reproduce the observed properties of 
the prompt GRB emission (Giannios & Spruit 2007; Giannios 2008). High efficiency of the 
conversion of the magnetic energy into the energy of radiating particles makes this model a 
viable alternative to the internal shock model. 

This work was supported by the US-Israeli Binational Science Foundation under grant 
number 2006170 and by the Israeli Science Foundation under grant number 737/07. 
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Fig. 1. — The development of the Kruskal-Schwarzschild instabihty in accelerating flows. 
In the presence of an effective gravity force, the plasma drips out of the current sheet that 
separates the oppositely directed magnetic fields. This promotes the magnetic reconnection. 



